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' I review the lattice formulations of vector-like gauge theories (e.g. QCD) with 

f*^ domain-wall fermions, and discuss how to optimize the chiral symmetry for any 

, finite N a (in the fifth dimension) , as well as to eliminate its dependence on as . 

o 

_f2 ' 1< Introduction 

Qh! A viable approach to study strongly-coupled gauge theories (e.g. QCD) is to 

formulate these theories on a spacetime lattice with domain-wall fermions. 



The basic idea of domain- wall fermions (DWF^P^is to use an infinite set of 
coupled Dirac fermion fields {ip s (x),s S (—00,00)} with masses behaving 
like a step function m(s) = m9(s) such that Weyl fermion states can arise 



> 
• 1— I 

x. 

as zeromodes bound to the mass defect at s = 0. However, if one uses a 
compact set of masses, then the boundary conditions of the mass (step) 
function must lead to the occurrence of both left-handed and right-handed 
chiral fermion fields, i.e., a vector- like theory. For lattice QCD with DWF^, 
in practice, one can only use a finite number (N s ) of lattice Dirac fermion 
fields to set up the domain wall, thus the chiral symmetry of the quark 
fields (in the massless limit) is broken. Obviously, the discretization in 
the fifth dimension also introduces the lattice spacing 05 into the theory. 
Presumably, only in the limit N s — > 00 and 05 — > 0, the correct effective 
4D theory with exact chiral symmetry can be recovered. Now the relevant 
question is how to minimize its dependence on iV s and 05. Since, in general, 
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they are independent parameters, it may happen that even in the limit 
N s — > oo, the massless quark propagator is exactly chirally symmetric but 
still has a strong dependence on a$ 0) . If this is the case, then the limit 
(05 — > 0) has to be taken before one can measure physical observables at 
any finite lattice spacing a, which of course, is difficult for any practical 
computations. 

It turns out that the conventional DWF action with open boundary 
conditions^ suffers from: (i) The chiral symmetry of the quark propagator 
is not optimal for any finite N s ; (ii) The quark determinant and propagator 
are sensitive to the lattice spacing 05, e ven for very large N s . 

In this talk, I present a formulatior02l of lattice QCD with DWF, in 
which the quark propagator is a5-invariant and has optimal chiral symmetry 
for any N s and background gauge field. 



2. Problems of the conventional domain-wall fermions 

First, we outline the basic features of DWF on the lattice. In general, given 
DWF action Af[ip,^j] with fermion fields {ip(x, s),^(x, s); s = 1,---,N S }, 
one can construct the quark fields q(x),q(x) from the boundary modes, and 
obtain the quark propagator 3, in a background gauge field as 

{q{x)q[y)) = jm&F*'' — = ( q)x - v {) 

where m q is the bare quark mass, 

_ 1+75S q _ n£i(i+°5g fl )-n£i(i-°5g 8 ) f9 s 
c i-^' n£i(i+**.)+n^(i 

and {H s , s = I, ■ ■ ■ , N s } are Hermitian operators which depend on Aj. 
For the conventional DWF, H s is the same for all s, and is equal to 

H = H w {2 + j 5 ar H w )~\ H w = l5 D w , (3) 

where D w is the standard Wilson-Dirac operator plus a negative parameter 
—mo (0 < mo < 2). In the limit N s — > 00, S becomes sgn(-ff), 

Um 5=^^ = ^, (4) 

and the quark propagator possesses exact chiral symmetry 

lim [(D c + m 9 )" 1 75 + J5{D C + m,)" 1 ] = . (5) 



a Here the color and Dirac indices are suppressed. 
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However, for the conventional DWF, the quark propagator still depends 
on 05 through H © . It turns out that the effects of a§ cannot be neglected 
even for very large N s . This is the essential difficulty encountered in lattice 
QCD calculations with the conventional DWF. So the relevant problem for 
lattice QCD with DWF is how to construct a DWF action such that its S 
operator is independent of 05, for any a, N s , and gauge background. 

Another difficulty of the conventional DWF is that it does not preserve 
the chiral symmetry optimally for any finite NM In other words, its S 
operator 

(1 + a 5 H) N > - (1 - a 5 H) N - 
(l + a 5 H) N ' + {l~a 5 H) ] 

is not the optimal rational approximation of sgn(iJ). The deviation of 
S(a§H) from sgn(if) can be measured in terms of 

yt{sgn(tf) - S(a 5 H)}Y 



S ^ = I!, 1J "1 - a 5 HR(4H 2 ) (6) 



o~(S) = max 

VY^O 



< max|sgn(?7) - S{r})\ , (7) 



rty 

where {77} are eigenvalues of a§H . Using the simple identity 

|sgn(x) - S(x)\ = |1 - Vx^R(x 2 )\ , S(x) = xR{x 2 ) (8) 
which holds for any x 7^ and S(x) = xR(x 2 ), we can rewrite ||7J| as 



a(S) < max 1 - y/rfRtf) (9) 
{if} 

where {n 2 } are eigenvalues of a\H 2 . For the conventional DWF, the r.h.s. 
of @ is not the minimum for any given N s , i.e., R{x 2 ) is not the optimal 
rational approximation of (x 2 ) -1 / 2 . Obviously, the problem of finding the 
optimal rational approximation S op t{x) = xR opt (x 2 ) of sgn(x) with x S 
[—Xmax, —Xmin] U [x m i n 7 x max ] is equivalent to finding the optimal rational 
approximation R opt (x 2 ) of (x 2 ) -1 / 2 with x 2 S [i^.^i^J. 

According to de la Vallee-Poussin's theorem and Chebycheff 's theorem^, 
the necessary and sufficient condition for an irreducible rational polynomial 

(n.m), \ p n X n + Pn-lX^ 1 -\ h p , . ^ „x 

g m i m + q m -ix m H h go 

to be the optimal rational polynomial of the inverse square root function 
x" 1 / 2 , < x m i n <x< Xmax is that <5(x) = 1 — r("' m ) (x) has n + m+ 2 
alternate change of sign in the interval [x m im %max]i an d attains its maxima 
and minima (all with equal magnitude), say, 

5{x) = -A,+A,---,(-l)" +m+2 A 
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at consecutive points (xi, i = 1, • • • , n + m + 2) 

In other words, if A n,m ^ satisfies the above condition, then its error 
a(r (n ^)= max 1 - y/E r^ n ' m \x) 

is the minimum amongall irreducible rational polynomials of degree (n, to). 

It is easy to show^l that R(x 2 ) © of the conventional DWF is not 
the optimal rational approximation for (cc 2 ) -1 / 2 . The optimal rational 
approximation for the inverse square root function was first obtained by 
Zolotare\0 in 1877, using Jacobian elliptic functions. A detailed discussion 
of Zolotarev's result can be found in Akhiezer's two bookg^l 

Thus the relevant problem for lattice QCD with DWF is to construct a 
DWF action such that the operator S in the quark propagator (JJJ is equal 
to 



, o f^iWK). N s = 2n + l, 

s opt - s opt (H w ) - <; HwR{ n-x, n){Hl) Ns = 2jh < 



where Rz(H^) is the Zolotarev optimal rational polynomial^! 



R 



z 



^{Hw) ~ > TT i . , 2 W / ' ~ Hw/^rnin (H) 

Amin h±/C->l-l 



1=1 



and 



p(w-l,w)/TT-2\ _ ^0 IL = 1 (1 + ^w/ C 2l) :,.), 



and the coefficients do , dp, ci and are expressed in terms of elliptic 
function|3 w ith arguments depending only on N s and b = tfnax/^min (Amin 
and Xmax are the minimum and the maximum of the eigenvalues of |-ff„,|). 



3. The optimal domain-wall fermions 

Recently, I have constructed a new lattice DWF actiorP^ such that the 
quark propagator is as-invariant and preserves the chiral symmetry op- 
timally for any N s and background gauge field. Further, its effective 
4D lattice Dirac operator for the internal fermion loops is shown to be 
exponentially-local for sufficiently smooth gauge background^. 
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Explicitly, the optimal lattice domain-wall fermion action reads b 

N s + 1 

A f = ^ ^ j>(x, s){(uj s a 5 D w (x, x') + 6 X<X >)5 S>S > 

s,s'—0 x,x' 

-{5*,x> ~ uj s a 5 D w (x,x'))(P-S s/ , s+1 +P+6 s , tS - 1 )}Tp(x',s') (13) 
with boundary conditions 

P + i/j(x,-1) = -rm q P + ij){x,N s + l) , r = (14) 

P_^(x, N s + 2) = -rm q P- ip(x, 0) , (15) 

and the quark fields constructed from the left and right boundary modes 

q(x) = [P-^{x, 0) + P+^x, N s + 1)] , (16) 
q(x) = s/¥ $(x 1 0)P + +4>(x,N s + l)P-] , (17) 

where the weights {lo s } are given by 

= un s +i = , (18) 
io s a 5 = y^v 7 ! - «' 2 sn 2 (v s ; k'), s = 1,---,N s . (19) 

Here sn(t! s ; k') is the Jacobian elliptic function with argument v s (|20() and 
modulus k' = y/1 - 1/6, where & = \ 2 max /\ 2 min . 

For the optimal DWF, H s = oj s H w in the operator 5 and the 
weights in (|19|l are obtained from the roots (u s = (w^as) -2 , s = 1, • • • , N s ) 
of the equation 



6z(u) = 



1 - ^uR { z' n) (u) = , iV s = 2n + 1 
1 - Vtti4 n ~ 1 '" ) ( u ) = . N s = 2n 



such that 5 is equal to S^pt pOfl- the optimal rational approximation of 
sga(H w ), and the quark propagator JQ) has the optimal chiral symmetry 
for any N s and b = \ 2 max /\ 2 mm . 
The argument v s in l|19|) is 

-=<- 1 »'^'^ 1 (\/M :VT ^) + [Sf (20) 



b In this paper, we suppress the lattice spacing a, as well as the Dirac and color indices, 
which can be restored easily. 
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where 

A=n — — * , m=tt — — ^- , (21) 

if' is the complete elliptic integral of the first kind with modulus k', and 9 is 
the elliptic theta function. From (|T§|> . it is clear that X^ ax < to s a 5 < A~ in , 
since sn 2 (; ) < 1. 

Obviously, the optimal domain-wall fermion action is invariant for 
any as, since its dependence on as is only through the product uj s a^ (|19|) 
which only depends on N s , A m ;„ and tfnax/^min- Therefore, 05 is a redun- 
dant parameter in the optimal domain-wall fermion action (|13|) . 

The generating functional for n-point Green's function of quark fields q 
and q is defined as 

[ ' 1 f[dU][d<ip][diJ}][d<i>}[dct)]e- A s- A f- A pf 1 ' 

where A ff is the gauge field action, Af is the domain- wall fermion action, 
A p f is the corresponding pseudofermion action with m q — 2mo, and J 
and J are the Grassman sources of q and q respectively The purpose 
of introducing the pseudofermion fields (which carry all attributes of the 
fermion fields but obey the Bose statistics) is to provide the denominator 
(1 + rDc)^ 1 in the effective 4D lattice Dirac operator 

Dim^^iDc + m^l + rDc)- 1 , D c = 2m ( ] + 75 f ° pt ) (23) 



\ 1 — "t$S op t j 

for internal fermion loops such that D{m q ) is exponentially-locaPI (for any 
N s and m q ), and its fermion determinant ratio (in the limit N s — > oo and 
a — > 0) is equal to the corresponding ratio of det[7 M (<9p + iA^) + m q \. 

Evaluating the integrations over {-0, -0} and {<fi, 0} in 122|) . one obtain:^ 

f{dU}e- A ^det[D(m q )}e J ( D c+™^J 
L ' J /[dE7]e- A «Mdst[I>(ro,)] 1 ' 

Then any n-point Green's function can be obtained by differentiating 
Z[J, J] with respect to J and J successively. In particular, the quark prop- 
agator is 

_ f{dU}e- A sdet[D(m q )}(D c + m q )-i y 
j=J=o j[dU]e- A ^et[D{m q )] 



S 2 Z[J,J] 



SJ(x)SJ(y) 
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which, in a background gauge field, becomes 

(q(x)q(y)) = (D c + m q )^ y , (25) 

and it goes to [jn{d n + iA^) + in the limit N s — + oo and a — > 0. 

For any gauge background and AT S , D c has the optimal chiral symmetry 
since the error a(S opt ) is the minimum, and a{S opt ) < (1 — A)/(l + A) ~ 
A(6)e-°( b ) iv % where A is defined in and c(6) can be estimated 

as A(b) ~ 4.06(l)6-°- 0091 ( 1 )ln(6)°- 0042 ( 3 ), c(6) ~ 4.27(45)ln(6)-°- 746 ( 5 ), b = 

\2 /\2 
^max / rain ' 

Finally, we note that D(0) (|23l) (N s — ► oo) times a factor r = l/2mo is 
exactly equal to the overlap Dirac operatoJ23l? = [1+D W (D\ JU D W )^ 1 / 2 ]I2. 
This implies that D is topologically-proper (i.e., with the correct index and 
axial anomaly), similar to the case of overlap Dirac operator. For any finite 
N s , r times D(0) (|23() is exactly equal to the overlap Dirac operator with 
(if 2 ) -1 / 2 approximated by Zolotarev optimal rational polynomial, (|1H (N s 
= odd) , or fT2^l ( N s = even) . Further discussions can be found in Ref. [6] . 
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